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Motivated by a recent experiment reporting on the possible application of graphene as sensors, 
we calculate transport properties of 2D graphene monolayers in the presence of adsorbed molecules. 
We find that the adsorbed molecules, acting as compensators that partially neutralize the random 
charged impurity centers in the substrate, enhance the graphene mobility without much change 
in the carrier density. We predict that subsequent field-effect measurements should preserve this 
higher mobility for both electrons and holes, but with a voltage induced electron-hole asymmetry 
that depends on whether the adsorbed molecule was an electron or hole donor in the compensation 
process. We also calculate the low density magnetoresistance and find good quantitative agreement 
with experimental results. 
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I. INTRODUCTION 

The recent discovery of graphene - a single layer of 
graphite - followed by the rapid progress in fabricating 
transistor-like devices and measuring its transport prop- 
erties^ is an important, perhaps seminal, development 
in low dimensional electronic phenomena in nanostruc- 
tures. These systems are conceptually novel, where the 
low-energy description for a single sheet of Carbon atoms 
in a honeycomb lattice is the linear "Dirac-like" disper- 
sion having both electron and hole carriers (which are the 
positive and negative chiral solutions of the Dirac Hamil- 
tonian) . The band structure induced carrier spectrum in 
graphene monolayers is four-fold degenerate (spin and 
valley), and the intrinsic sublattice symmetry (two in- 
equivalent Carbon atoms in the unit cell) causes a sup- 
pression of backscattering. Both of these features could 
have application in technology where the former pro- 
vides the opportunity to have both spin-tronic and valley- 
tronic functionality on the same device, and the theoret- 
ical absence of backscattering has led to the speculation 
that, as a matter of principle, carrier mobilities in 2D 
graphene monolayers could be extremely high even at 
room temperature. 

The increase in graphene carrier mobility induced by 
absorbed gas molecules has been recently used as a 
highly-sensitive solid-state sensor capable of detecting in- 
dividual molecules^ Notwithstanding the interpretation 
made in Ref. the experimental data reported there 
certainly shows Hall mobility /i# — Pxyf(PxxB) increas- 
ing with dopant exposure time. The experimental data 
is reproduced in Fig. [1] where for both NH3 and NO2 
adsorbants one finds a strong increase in mobility fin for 
small doping concentrations and that this increase in mo- 
bility persists as one increases the dopant exposure time 
t. It is further reported in Refs.[i][5|that the hall mobility 
then saturates for even larger doping concentrations or 
during subsequent annealing of the samples. 

While this experimental advance has the potential to 
revolutionize gas sensors, it also raises fundamental ques- 
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FIG. 1: (Color online) Hall mobility \u,h | in units of m 2 /Vs 
measured in the experiments of Ref. |3| by chemically doping 
graphene with NH3 and NO2. Notice that \[ih\ increases as 
one increases the dopant exposure time t, indicating a de- 
crease in impurity scattering. 



tions about the mechanism by which adsorbates change 
the transport properties of graphene. Similar exper- 
iments done on Carbon nanotubes&i interpreted their 
data to argue that NH3 adsorbents (in the presence of 
water) transfered 0.04e~ per molecule, while NO2 binds 
to the surface and withdraws 0.1e~ per molecule. 8 How- 
ever, these experiments on Carbon Nanotubes had no di- 
rect way to measure the carrier density, and were unable 
to provide further clues as to the charge transfer mecha- 
nism. In contrast, the recent experiments of Ref. [HIH on 
"chemically-doped" graphene (which is a zero-gap semi- 
conductor with both electron and hole carriers) could 
provide definitive answers to these long standing theo- 
retical questions. Although there has been considerable 
recent theoretical activity^ ' 10 ! 11 ' 12 ! 13 ' 14 ! 15 studying car- 
rier transport in graphene, transport in chemically doped 
graphene - the subject of our current work - has not been 
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FIG. 2: (Color online) Graphene magnetoresistivity p xx and 
p xy for B — IT. Dashed (solid) lines are calculated for one 
(two) carrier model, Eq. [T]with impurity concentration m = 
1.75xl0 12 cm -2 . Open circles show data of Ref. OJH Inset 
compares longitudinal conductivity a xx for the two carrier 
model (solid line) and experiments of Ref. 00 (circles). 



considered previously in the literature. Since it is now 
believe d 16 ! 17 that charged impurities are the dominant 
scattering mechanism in graphene, one of the main goals 
of the present work is to understand chemical doping of 
graphene within the framework of charged impurity scat- 
tering. 

We emphasize that the carrier mobility n is the average 
velocity of each electron per unit applied electron field. 
It is strictly defined only in the limit of vanishing elec- 
tric and magnetic fields. However, since for small fields, 
the electrical conductivity is the ratio of the current J to 
the electric field E, we have the relation between conduc- 
tivity and mobility as a = neu. Therefore, the correct 
way to obtain the mobility from the measured conductiv- 
ity is to divide the conductivity by a carrier density, i.e. 
fi = cr/ne. Only in the very special case where conduc- 
tivity varies strictly linearly with carrier density is this 
result the same as the derivative of conductivity with 
respect to density, i.e. (l/e)da/dn, but in general the 
mobility can be quite different from (l/e)da/dn when 
the conductivity is not a linear function of density. To 
our knowledge, the theoretical justification for a linear- 
in-density conductivity is valid only for incoherent trans- 
port of uncorrelated Coulomb impurities located exactly 
at the surface of the graphene sheet. Other scattering 
mechanisms likely to be present in graphene including 
phonons, ripples, and defects are known to change the 
density dependence of the conductivity necessitating the 
consistent use of fj, = cr/ne when discussing carrier mo- 
bility and especially when discussing multiple scattering 
mechanisms. This is relevant in comparing our theory 
with the data of Ref. where the mobility has been de- 
fined as (l/e)da/dn rather than the conventional a/ne. 

We note that for most graphene samples, the mea- 



sured conductivities are linear over wide range of density. 
This fact alone is the smoking gun for the dominance 
of charged impurity scattering in current graphene sam- 
ples, although we note that other signatures of charged 
impurity scattering have been observed in recent exper- 
iments i 16 ! 17 We note that even within the charged im- 
purity scattering model, the linear-in-density behavior 
breaks down at low densities near the Dirac point, where 
instead one expects a constant conductivity plateau 
whose value is determined by the residual density in- 
duced by the charged impurities*^ At high densities (i.e. 
high gate voltages) one expects^ the competing effects 
of short-range scattering (such as point-defects or dislo- 
cations) and shifting of the charged impurities that also 
cause the conductivity to deviate from the the linear re- 
lation. For all these reasons, to get the correct mobility 
one has to divide the conductivity by the density. 

In a finite magnetic field, the Hall mobility is often used 
to characterize samples and is obtained by measuring the 
conductivity a and then dividing it by the concentration 
n found by the Hall effect. The Hall mobility has the 
same dimension as the drift mobility discussed above, 
but differs by a factor ru , which is the scattering factor 
of the Hall coefficient, defined by rn = ( t2 )/{ t ) 2 , where 
t is a relaxation time and (• • • ) indicates the energy aver- 
age. For a model with only charged impurity scattering 
and away from the Dirac point and at zero temperature, 
one can calculate that the difference between the defini- 
tions of Hall mobility and drift mobility is less than 3 
percent, making one a good approximation for the other. 
But again, this is only true away from the Dirac point, 
and in the case of charged impurity scattering being the 
dominant scattering mechanism. In general, to get the 
correct Hall mobility experimentally, one would have to 
divide the measured Hall conductivity by the carrier den- 
sity instead of doing a derivative. 

We now proceed as follows. In Section [H] we first gen- 
eralize the theory^ of charged impurity scattering to in- 
clude the effect of a magnetic field and compare our the- 
ory for graphene magnetoresistance with the experimen- 
tal results of Ref. d in the absence of chemical dopants. 
Having understood the zero doping regime, in Section ITVl 
we discuss the most natural assumption for the effect of 
chemical adsorbents which is to transfer charge to the 
graphene layer leaving behind a charged impurity. We 
find that this always decreases the graphene mobility in 
contrast to the experimental findings in Ref. 3j. In Sec- 
tion |IIT] we propose a compensation model to explain the 
increase in mobility seen experimentally. In addition, 
we highlight that this model makes several predictions 
that can be verified experimentally. In Section [V] we 
phenomenologically account for the observed^ saturation 
of mobility during annealing experiments on molecular 
doped graphene sheets. We note that in our model a 
saturation of mobility necessarily means that there is 
no change in the charged impurity concentration and 
we demonstrate that including an asymmetric shift in 
threshold voltage reproduces all the features observed 
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experimentally, 
a discussion. 



We then conclude in Section I VII with 



II. GRAPHENE MAGNETORESISTANCE 

Before considering chemical doping issues, we first dis- 
cuss graphene magnetoresistance. The transport proper- 
ties of graphene in the presence of a magnetic held was 
considered for "white noise" (i.e. delta-correlated) dis- 
order by Zheng and Ando^, and in the Landau quan- 
tization limit by Gusynin and Sharapovi 2 ^ As has been 
demonstrated theoretically in Refs. [9l lTolfl3l . and exper- 
imentally in Refs. flMl7l . the dominant transport mech- 
anism in graphene is Coulomb scattering from charged 
impurities (not white noise disorder), and below we ex- 
tend this formalism to incorporate the effects of a finite 
(but weak) magnetic field using the Boltzmann transport 
theory in the semiclassical regime where Landau quan- 
tization is unimportant. While we are concerned with 
the high density limit, where n > n, and where Boltz- 
mann theory is valid, to compare with experiments we 
construct a simplified two-component (i.e. electron and 
hole) model consistent with the percolation model sug- 
gested in Ref. U To phenomenologically account for the 
conductivity saturation at low density within the Boltz- 
mann theory framework, we assume that for V g > 0, the 
electron and hole densities are 



n e — tiq + n, 

n h = {n Q - n)0(n o - n), 



(1) 



respectively, where 6{x) is the Heavyside step function 
and no ~ rii is the minimum carrier density. We have 
similar expressions for V g < 0. In this simplified model, 
uq — for the one carrier model and we use 2no = n,-/4 
for the two carrier models Here, we calculate the mo- 
bility in the presence of randomly distributed Coulomb 
impurity charges near the surface with the electron- 
impurity interaction being screened by the 2D electron 
gas in the random phase approximation (RPA). The 
screened Coulomb scattering is the only important scat- 
tering mechanism in our calculation. We assume that the 
direction of current flow is in the x direction and that the 
magnetic field is in the z direction, with the graphene 
layer being in the 2D xy plane. In the presence of two 
types of carriers (electrons and holes) in a finite magnetic 
field B, the current density in the x and y directions are 
given by 



Jy = {o$+v$]E x + [a$+a£J]E y . 
The longitudinal conductivities are given by 



r(<0 



a^R^B 



(2) 
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where the superscript c = e,h denotes electron and hole 
carriers, and the Hall conductivities are given by 
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Here the zero-field electrical conductivity for each 



carrier is defined by a, 



2 'MIE^t/K) 



20(e 2 /h)(n {c) /ni) (see Refs. iQMlllll for details), 



and the Hall coefficient 



l/« (c) e (c) , 



where rir c -\ 



n e ,rih and — ±e, are the density and charge of the 
carriers, respectively. 

Shown in Fig. [2] is the Boltzmann magnetoresistivity 
theory (Eqs.OS]) for the two carrier model that is com- 
pared with the gate-voltage tuned experimental results 
reported in Ref. d without any chemical doping. The 
remarkable agreement with experimental data is further 
highlighted in the inset where the theoretical longitudinal 
conductivities are compared. There is no free parameter 
in the theory (although experimental results have been 
scaled by rii, determined from the conductivity depen- 
dence at high density). For comparison, dashed lines in 
the main panel of Fig. [2] show resistivities with one kind 
of carrier. Since in general, the conductivities (a xx and 
<J xy ) are proportional to the carrier density n, the calcu- 
lated resistivities for the single component model diverge 
as n — > 0. However, in two component model, the Hall 
conductivity is given by a xy = a xy + a xy and a xy oc R l H . 
Since the Hall coefficient R l H has different signs due to the 
sign of charge, the sum of the Hall conductivity becomes 
zero when n e — nh- This results in a finite longitudinal 
magnetoresistivity and a zero Hall resistivity at n = 0. 
Similar to the case of zero magnetic field, for finite mag- 
netic field, the high density mobility is the only reliable 
measure of the concentration of Coulomb scattering cen- 
ters. For graphene on a SiC>2 substrate we find 



'xy 



p xx B <J XX B 

48.37 f n e - n h \ 48.37 

5 + n h 



(5) 



where impurity density ni is measured in units of 
10 10 cm -2 and hh is measured in units of m 2 /Vs. 



III. DECREASED MOBILITY REGIME 

We first consider the most natural assumption for the 
effect of chemical adsorbents which is to transfer charge 
to the graphene layer leaving behind a charged impurity 
at some distance d ~ 5 A from the surface. This is the 
"standard model" assumed, for example, in carbon nan- 
otubes. 8 Shown in Fig. [3] are the results of a Boltzmann 
calculation for this case. The calculation shown in the 
left panel assumes "conservation of charge" , i.e. for ev- 
ery carrier induced by the adsorbent there is left behind 
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FIG. 3: (Color online) Left panel shows calculated density 
dependent resistivity for different locations of the impurity 
layer d assuming that adsorbed chemicals contribute to the 
carrier density n leaving behind a charged impurity at 
a distance d from the graphene surface. Right panel shows 
calculated resistivity using the same model (with d — 5 A) 
but assuming that only a fraction of the transferred charge 
contributes to the impurity density. 



an impurity of equal charge that on the average is at a 
distance d from the graphene surface. On the right panel 
we relax this assumption and for fixed d = 5 A, we show 
calculated results for the impurity charge tin being n/2 
or n/3, where n is the induced carrier density. This model 
would always cause: (i) a decrease in mobility, because 
of the increased impurity scattering, and (ii) an increase 
in carrier density, because of the charge transferred by 
the adsorbed molecule. However, both these conclusions 
are at odds with recent experiments^ This forces us to 
conclude that the dominant source of scattering remains 
the random charged impurities at the surface between 
graphene and the S1O2 substrate, and that the chemical 
adsorbent compensates some of this charged background 
thereby providing an increased mobility. 



IV. INCREASED MOBILITY REGIME 

Within our magnetoresistance model, an increase in 
Hall mobility (for large carrier density) necessarily im- 
plies a decrease in the charged impurity concentration, 
which we argue is a consequence of adsorbed molecules 
acting as compensators that partially neutralize the ran- 
dom Coulomb scatterers in the substrate. We note that 
molecular dopants may also shift the threshold voltage as 
a consequence of surface chemistry effects and this pro- 
vides no information about the density of charged impu- 
rities contributing to graphene transport. An important 
goal of the current work is to theoretically explore the 
consequences of the compensation assumption using the 
full Boltzmann transport theory and to make predictions 
for subsequent field-effect conductivity measurements on 



the chemically doped graphene sheets. Since the prospect 
for applications of graphene as sensors is an important 
technological possibility, our work has practical implica- 
tions. More importantly, our work would critically vali- 
date (or invalidate) the transport model and mechanism 
in graphene 

From the data reported in Ref.d, one can conclude that 
initial chemical doping serves to neutralize the charged 
impurities in the substrate while not contributing to the 
carrier density. The experimental data motivates the as- 
sumption that there is a roughly equal number of posi- 
tive and negative charged impurities so that nf ~ . 
This mechanism of adsorbed molecules acting as compen- 
sators should eventually saturate, after which the mobil- 
ity would decrease as shown in Fig. [3l Before doping, the 
interface between graphene and the substrate has a fixed 
density of random charged impurities rij = nf , com- 
prising both positively (nf) and negatively (n~) charged 
impurities. The electron-hole asymmetry seen in the ex- 
periments could be caused by the positive (negative) gate 
voltage shifting positively (negatively) charged impuri- 
ties by a distance d away from the surface. We take this 
effect into account explicitly, but unlike Ref. [^, we as- 
sume for simplicity that before doping nf — n~ . We 
also assume that nf = for graphene doped with NH 3 
(NO%). This model provides definitive predictions that 
could be tested in future experiments. For example, in 
the majority of samples where point (i.e. short-range) 
scattering is unimportant, chemical doping should in- 
crease the electron- hole asymmetry with NO2 (NH3) 
showing conductivity as a function of carrier density for 
electrons (holes) that increases faster than the linear be- 
havior predicted for holes (electrons). In the presence of 
both charge and point scattering, applying a gate voltage 
results in a non-universal crossover between linear and 
constant conductivity. This correction caused by short- 
range scatterers has been called the "sub-linear" con- 
ductivity, in contrast to the "super-linear" conductivity 
described above to characterize the electron-hole asym- 
metry. For samples with sub-linear behavior, chemical- 
doping would shift the cross-over closer to the constant 
conductivity which would be seen experimentally as shift 
of the onset of the sub-linear conductivity to lower den- 
sity. If these features are seen in future experiments on 
graphene, they would rule out alternative explanations 
including the conventional wisdom in the nanotube com- 
munity that chemical adsorbents contribute directly to 
the graphene carrier density rather than neutralize the 
substrate^, as we claim here. 

The predictions of our theoretical model are shown in 
Figs. |U First, we predict that because chemical dop- 
ing neutralizes the interface impurities, the high mobil- 
ity should persist when doing field-effect measurements 
after exposure to the adsorbant gases. This can be seen 
in Fig.[4l where the chemically-doped cases have a higher 
slope than the undoped case. Assuming that point- 
scattering is still unimportant, the slope should continue 
to be linear as in the undoped case. Second, we predict a 
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FIG. 4: (Color online) Shown (solid lines) are the calculated 
graphene conductivity for both doped and undoped cases for 
m = 1.75 x 10 12 cm" 2 . The asymmetry between NH3 doping 
and NO2 doping is caused by the voltage induced shifting of 
impurities as explained in the text. Also shown (dashed line) 
is the undoped case without any voltage induced impurity 
shifting. Note that for both NH3 and NO2 , doping increases 
the sample mobility. 

slight asymmetry between hole doping (NO2) and elec- 
tron doping (NH3) that is caused by the same mechanism 
described above to explain the electron-hole asymmetry 
in the undoped case. The physical picture underlying our 
transport model for graphene is simple: carrier transport 
in 2D graphene layers is determined by charge impurity 
scattering with the possibility of partially suppressing the 
impurity scattering through adsorbate-induced compen- 
sation. At low carrier density, both electrons and holes 
are present, leading to the observed^ magnetoresistive be- 
havior. If short-range point scatterers are present in the 
system in addition to charged Coulomb scatterers then 
the mobility will eventually "saturate" (i.e. become sub- 
linear in carrier density) at high enough carrier densities. 

V. SATURATION OF MOBILITY 

We observe that our two-carrier model also accounts 
phenomenologically for the recent annealing experiments 
on molecular doped graphene sheets reported in Ref. la- 
in contrast to Ref. [H, those experiments show no change 
in the mobility implying that there is no change in the 
charged impurity concentration n.j. The experiments 
however do show a shift in threshold voltage, which 
can be easily modeled by replacing n in Eq. Q] with 
h = n — (1 =F v)^ f° r electrons (holes). Using 77 w 0.2 to 
characterize the asymmetry in the data between electrons 
and holes, we find that this reproduces all the features 
of the experiment including that for electrons, the l/p xy 
traces are shifted in parallel to lower values, that the peak 
corresponding to the divergence in l/p xy is also shifted 
to higher gate voltages, and is broader and its minimum 
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FIG. 5: (Color online) Main panel shows two carrier model 
used to explain the results of Ref. 0. Annealing causes a shift 
in the threshold voltage parameterized by (from left to right) 
A/n = 0,3,6,9. 

value increased £ These features are shown in Fig. We 
point out that such features involving different devices 
having different threshold voltages but identical mobili- 
ties at high densities are quite common in Si MOSFETs 
and GaAs 2D hcterostructures. 



VI. CONCLUSION 

We have generalized the charged impurity scattering 
model (Ref. M) of graphene transport to calculate the 
low-field magnetoresistance and the chemical doping sit- 
uation as explored in recent experiments We obtain 
reasonable agreement with experiments by assuming that 
the chemical dopants lead to two distinct effects: increas- 
ing carrier density and decreasing (i.e. compensating) the 
already existing charged impurity density. The mobility 
is a function of three parameters in the system: carrier 
density (n), impurity density (n^), and impurity distri- 
bution (which can be parameterized in a simple model 
by a single length scale d giving the separation of the 
charged impurities in the substrate from the graphene 
carriers. In our model, the effect of chemical doping on 
mobility must enter through the modification of these 
three lengthscales: c?, n -1 / 2 , n" 1 ^ 2 . How chemical dop- 
ing affects these three parameters microscopically is not 
known for the experiments of Refs. We have there- 
fore made the simplest possible assumption of keeping 
the unknown separation parameter d a constant, assum- 
ing that chemical doping can only affect n and n^. In 
principle, it is possible for chemical doping to modify 
the impurity separation d, giving us another handle on 
the modification of mobility by chemical doping. The 
important point to emphasize here is that to the extent 
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chemical doping increases (decreases) the charged impu- 
rity density rii, it must change mobility/i ~ n' 1 . If /i 
remains unaffected as claimed in Refs.|j|a, then chemical 
doping must somehow affect only n without affecting d 
and rii (the other possibility is a chemical doping induced 
complex interplay among n, rij, d and perhaps even the 
background dielectric constant k, conspiring to keep mo- 
bility a constant, a situation clearly beyond the scope 
of any minimal theory). We believe that the theory pre- 
sented herein provides the most straightforward explana- 
tion of the experimental observations within the minimal 
model (and keeping the number of adjustable parame- 
ters to a minimum). It is gratifying in this context to 
point out that a recent experimental studyi£ of system- 
atic potassium doping effects on graphene mobility in 
ultra-high vacuum obtains excellent agreement with our 
minimal theory-SilS, using only m as the single unknown 
parameter. 

In summary, we have developed a detailed micro- 



scopic theory for graphene carrier transport which ex- 
plains very well the existing magnetoresistance data as 
well as accounting qualitatively and quantitatively for 
the recently observed adsorbate-induced modification of 
graphene transport properties which has implications for 
application as sensors. We make a number of specific 
experimental predictions based on our transport model, 
whose validation (or falsification) should further con- 
solidate our understanding of 2D carrier transport in 
graphene. Our theory provides quantitative estimates for 
the constraints on using graphene transport properties as 
chemical sensors or detectors. 
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